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Beschreibung 

Retrospektive bzw. f enstergesteuerte Filterung von Bildern 
zur Adaption von Scharfe und Rauschen in der Computer- 
5 Tomographic 

Mit modernen medizinischen Diagnoseverf ahren, wie beispiels- 
weise der Rontgen-Computertomographie (CT) , konnen Bilddaten 
eines untersuchten Messobjektes gewonnen werden. In der Regel 
10 handelt es sich bei dem untersuchten Messobjekt urn einen 
Patienten . 

Die Rontgen-Computertomographie - im Nachf olgenden kurz mit 
CT bezeichnet - ist ein spezielles Rontgen-Aufnahmeverf ahren, 

15 dass sich im Bildaufbau grundsatzlich von dem klassischen 

Rontgen-Schichtauf nahmeverf ahren unterscheidet . Bei CT-Auf- 
nahmen erhalt man Transversal-Schnittbilder , also Abbildungen 
von Korperschichten, die im Wesentlichen senkrecht zur Kor- 
perachse orientiert sind. Die im Bild dargestellte gewebespe- 

20 zifische physikalische Grofle ist die Verteilung des Schwa- 

chunswertes von Rontgenstrahlung |X(x,y) in der Schnittebene . 
Das CT-Bild erhalt man durch Rekonstruktion der vom verwende- 
ten MeiJsystem gelieferten eindimensionalen Projektionen der 
zweidimensionalen Verteilung von |i(x,y) aus zahlreichen 
^5 verschiedenen Blickwinkeln . 

CT-Bilder sind sowohl mittels eines CT-Gerats mit endlos 
umlauf f ahigem Abtastsystem als auch mittels eines C-Bogen- 
Gerats mit nur urn weniger als 360° rotierbarem Abtastsystem 
30 generierbar. Die Abkurzung "CT" , z.B. bei "CT-Rohdaten" , wird 
im Folgenden im Hinblick auf beide Geratetypen gebraucht . 

Die Projektionsdaten ermitteln sich aus der Intensitat I 
eines Rontgenstrahls nach seinem Weg durch die abzubildende 
35 Schicht und seiner ursprunglichen Intensitat I 0 an der Ront- 
genquelle gemali dem Absorptionsgeset z 
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Der Integrationsweg L reprasentiert die Bahn des betrachteten 
Rontgenstrahls durch die zweidimensionale Schwachungsvertei- 
5 lung |i(x,y). Eine Bildpro j ektion setzt sich dann aus den mit 
den Rontgenstrahlen einer Blickrichtung gewonnenen Messwerten 
der Linienintegrale durch die Obj ektschicht zusammen. 

Man erhalt die aus verschiedensten Richtungen stammenden 
10 Projektionen - charakterisiert durch den Pro j ektionswinkel a 
- durch ein koinbiniertes Rontgenrohren-Detektor-System, das 
in der Schichtebene urn das Objekt rotiert. Die derzeit ge- 
brauchlichsten Gerate sind sogenannte „Facherstrahlgerate* 
bei denen Rohre und ein Array aus Detektoren (eine lineare 
15 Anordnung von Detektoren) in der Schichtebene gemeinsam um 
ein Drehzentrum, welches auch Mitte des kreisf ormigen Mess- 
feldes ist, rotieren. Die mit sehr langen Messzeiten behafte- 
ten „Parallelstrahlgerate* werden hier nicht erlautert. Es 
sei jedoch darauf hingewiesen, dass eine Transformation von 
20 Facher- auf Parallelpro jektionen und umgekehrt moglich ist, 
so dass die vorliegende Erfindung, die anhand eines Facher- 
strahlgerates erklart werden soli, ohne Einschrankung auch 
fur Parallelstrahlgerate anwendbar ist. 

^5 Bei Facherstrahlgeometrie besteht eine CT-Aufnahme aus Li- 

nienintegralmesswerten -ln(I/I 0 ) eintref fender Strahlen, die 
durch eine zweidimensionale Verkniipfung des Pro j ektionswin- 
kels a€ [0,27i) und den die Detektorpositionen def inierenden 
Facherwinkeln (ie[-f3o,Po] (Po ist der halbe Facherof f nungswin- 

30 kel) charakterisiert sind. Da das Messsystem nur liber eine 
endliche Anzahl k von Detektorelementen verfiigt und eine 
Messung aus einer endlichen Anzahl y von Projektionen be- 
steht, ist diese Verkniipfung diskret und kann durch eine 
Matrix dargestellt werden: 

35 



p(*y,Pk)- [0,271) x [-p 0 ,Po] 
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bzw . 



p (y,k): (1,2,... 



N P ) x (1,2 



, . . . 



N s ) 



(3) 



Die Matrix p CMO heiBt Sinugramm fur Facherstrahlgeometrie . 
Die Projektionszahl y und die Kanalzahl k liegen in der 
Groflenordnung von 1000. 

10 Bildet man gemafi Gleichung (1) die Logarithnien, so erhalt man 
also die Linienintegrale samtlicher Projektionen 



15 deren Gesamtheit auch als die Radon-Transf ormierte der Ver- 
teilung \x(x,y) bezeichnet wird. Eine solche Radon- 
Transformation ist umkehrbar, pi(x,y) demnach aus p(a,(3) durch 
Rucktrans formation (inverse Radontrans format ion) errechenbar . 
Bei dieser Rucktransf ormation wird ublicherweise ein Fal- 

20 tungs-Algorithmus angewendet, bei dem die Linienintegrale je 
Projektion zunachst mit einer speziellen Funktion gefaltet 
und dann langs der ursprunglichen Strahlrichtungen auf die 
Bildebene ruckpro j iziert werden. Diese spezielle Funktion, 

^ durch die der Faltungs-Algorithmus im Wesentlichen charakte- 

25 risiert ist, wird als " Faltungskern" bezeichnet. Durch die 

mathematische Gestalt des Fal tungskerns besteht die Moglich- 
keit, die Bildqualitat bei der Rekonstruktion eines CT-Bildes 
aus den CT-Rohdaten gezielt zu beeinf lussen . Beispielsweise 
konnen durch einen entsprechenden Faltungskern hohe Frequen- 

30 zen betont werden, urn die Ortsauf losung im Bild zu erhohen, 
oder durch einen entsprechend anders gearteten Faltungskern 
hohe Frequenzen zur Reduktion des Bildrauschens gedampft 
werden. Zusammengef asst kann also gesagt werden, dass bei der 
Bildrekonstruktion in der Computertomographie durch Wahl 

35 eines geeigneten Faltungskernes die Bildcharakteristik, die 
durch Bildscharf e/Bildrauschen und Bildkontrast (beide ver- 




/ 



(2) 
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halten sich zueinander komplementar ) charakterisiert ist, 
beeinflusst werden kann. Dabei besteht ein direkt proportio- 
naler Zusammenhang zwischen Bildscharfe und Bildrauschen, 
d.h. dass bei Erhohung der Bildscharfe im gleichen Maiie das 
5 Rauschen erhoht wird. 

Auf das Prinzip der Bildrekonstruktion in der CT durch Be- 
rechnung der |j,-Wert-Verteilung soil nun nicht weiter einge- 
gangen werden. Eine ausfiihrliche Beschreibung der CT- 
10 Bildrekonstruktion ist beispielsweise in „Bildgebende Systeme 
fur die medizinische Diagnostic , 3 .Auf lage, Munchen: Publicis 
MCD Verlag, 1995, Hrsg.: Morneburg Heinz, ISBN 3-89578-002-2, 
dargestellt . 

15 Allerdings ist mit der Berechnung der ja-Wert-Verteilung der 
durchstrahlten Schicht die Aufgabe der Bildrekonstruktion 
noch nicht abgeschlossen . Die Verteilung des Schwachungskoef- 
fizienten |i reprasentiert im medizinischen Anwendungsbereich 
nur eine anatomische Struktur, welche noch in der Form eines 

20 Rontgenbildes dargestellt werden mufl . 

Nach einem Vorschlag von G.N. Hounsfield ist es allgemein 
ublich geworden, die Werte des linearen Schwachungskoef f i- 
zienten \i (der die MaBeinheit cm" 1 hat) auf eine dimensions- 
^5 lose Skala zu transf ormieren, in der Wasser den Wert 0 und 
Luft den Wert -1000 erhalt. Die Umrechnungsf ormel auf diese 
„CT-Zahl TV lautet : 

CT-Zahl= M ~ MlVasser \000 (4) 
^Wasser 

30 

Die Einheit der CT-Zahl heifit „Hounsf ield-Unit^ (HU) . Diese 
als " Hounds f ield-Skala" bezeichnete Skala ist sehr gut zur 
Darstellung anatomischen Gewebes geeignet, da die Einheit HU 
die Abweichung in Promille von |a W a 3 ser ausdrlickt und die |a- 
35 Werte der meisten korpereigenen Substanzen sich nur wenig vom 
|j.-Wert des Wassers unterscheiden . Aus dem Zahlenbereich (von 



200201333 



5 

-1000 fur Luft bis ca 3000) werden nur ganze Zahlen als 
Trager der Bildinf ormation verwendet. 

Allerdings wiirde die Darstellung des gesamten Skalenbereiches 
5 von etwa 4000 Werten das Unterscheidungsvermogen des mensch- 
lichen Auges bei weitem ubersteigen. Zudem interessiert den 
Betrachter oft nur ein kleiner Ausschnitt des Schwachungswer- 
tebereiches, z.B. die Dif f erenzierung von grauer und weiiier 
Gehirnsubstanz, die sich nur urn etwa 10 HU unterscheiden. 

10 

Aus diesem Grunde benutzt man die sogenannte Bildf ensterung 
(engl.: Windowing). Dabei wird nur ein Teil der CT-Werteskala 
ausgewahlt und uber alle verfugbaren Graustufen gespreizt. 
Auch kleine Schwachungsunterschiede innerhalb des gewahlten 
15 Fensters werden so zu wahrnehmbaren Grautonunterschieden, 

wahrend alle CT-Werte unterhalb des Fensters schwarz und alle 
CT-Werte oberhalb des Fensters weift dargestellt sind. Das 
Bildfenster kann sowohl in seinem Zentralniveau als auch in 
seiner Weite beliebig variiert werden. 

20 

Nun ist es in der Computer-Tomographie von Interesse, Mul- 
tiplanar Reformatierungen (MPR bzw. Sekundarschnitte) vorzu- 
nehmen . Multiplanare Reformatierungen sind aus einem Volumen- 
datensatz (auch als Primardatensatz bezeichnet und 
^5 iiblicherweise reprasentiert durch diinne axiale Schichten) 

berechnete, beliebig geneigte, ebene CT-Bilder. Da die Pixel 
im Allgemeinen nicht die in dem Volumendatensat z definierte 
Position treffen und die Schichtdicke einer MPR beliebig 
eingestellt werden konnen soil, muss dabei geeignet interpo- 
30 liert werden. Die Auflosung eines Volumendatensatzes insbe- 

sondere bei neueren CT-Geraten ist nahezu isotrop. Aus diesem 
Grund lassen sich aus einem solchen Volumendatensat z hochqua- 
litative MPR berechnen deren Qualitat sich nicht von der der 
primaren Bilder unterscheidet . 

35 

Allerdings ist es im Rahmen einer guten diagnos tischen Bild- 
bewertung ebenfalls von Interesse, die Bildcharakteristik 
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eines MPR - im Wesentlichen gekennzeichnet durch Scharfe und 
Rauschen - durch geeignetes Filtern zu Manipulieren . Scharfe 
und Rauschen der Ref ormatierungen sind im Wesentlichen durch 
Scharfe und Rauschen der primaren, axialen Bilder sowie durch 
5 die bei der Generierung der MPR eingestellten Schichtdicke 
bestimmt . 

Die Darstellung von CT-Bildern, insbesondere die von MPRs, 
mit unterschiedlicher Bildcharakteristik ist deshalb von 
10 Interesse, weil unterschiedliches Evaluieren (d.h. eine 

unterschiedliche klinische Bewertung) der gleichen Aufnahme 
des entsprechenden Gewebes eine verschiedenartige Darstellung 
~ des auf genommenen Gewebes erfordert. 

15 Im Stand der Technik wird eine gezielte Manipulation der 

Bildcharakteristik von Sekundarschnitten dadurch erreicht, 
dass aus den Rohdaten ein neuer Volumendatensat z durch eine 
neue Bildrekonstruktion mit geanderten Faltungskernparametern 
ermittelt wird und auf der Basis dieses neuen Primardatensat- 

20 zes eine anschlieUende erneute Erzeugung der ursprlinglichen 
Sekundarschnitte erfolgt. Dies bedeutet eine Rekonstruktion 
mit einem anderen Faltungskern, der eine unterschiedliche 
Charakteristik aufweist, z.B. eben eine andere Scharfe. 

fj^> Durch diese bekannte Vorgehensweise kann das oben beschriebe- 
ne Ziel, die Bildcharakteristik von CT-Bildern insbesondere 
von Sekundarschnitten gemafr den diagnostischen Anf orderungen 
optimal anzupassen, nur zum Teil realisiert werden. Eine 
erneute erste Rekonstruktion mit den dem spateren Nutzen 
30 angepassten Parametern verandert die Scharfe nur transversal 
zur Patientenachse . Die Patientenachse wird im Folgenden als 
z-Achse bezeichnet. In einer aus dieser Rekonstruktion er- 
zeugten MPR bleiben die bildcharakterisierenden Eigenschaf ten 
in z-Richtung unverandert. 



Um auch in z-Richtung den gewiinschten Erfolg zu erzielen, 
mussten die Schichtprof ile (engl.: Slice Sensity Profile, 
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10 



15 



SSP) dieser ersten Rekonstruktion individuell angepasst 
werden. Bisher ist diese Anpassung nur in diskreten Schritten 
moglich, das heiftt eine Funktion zur beliebig genauen Abstim- 
mung der SSP - die im Prinzip auch als axiale Scharfe be- 
zeichnet werden kann - wird in derzeitigen CT-Geraten nicht 
zur Verfiigung gestellt. Auch wenn ein beliebiges SSP bei der 
ersten Rekonstruktion verwendet werden konnte, wurde die 
Erzeugung von geeignet veranderten Sekundarschnitten einen 
erheblichen Bedienungs- und Rechenzeitauf wand erfordern. 

Aufgabe der vorliegenden Erfindung ist es daher, hinsichtlich 
diagnostisch optimierten Gesichtspunkten neue Techniken zur 
Verbesserung bzw. Vereinf achung der Filterung von CT-Bildern 
vor zuschlagen . 

Erf indungsgemafi wird also ein erstes erf indungsgemaftes Ver- 
fahren zur Filterung von CT-Bildern vorgeschlagen, aufweisend 
die folgenden Schritte: 

20 a) Aufnehmen eines CT-Rohdatensatzes mit einem CT-Gerat oder 
mit einem C-Bogen-Gerat, 

b) Rekonstruktion eines Primardatensatzes aus dem CT- 
Rohdatensatz mittels einem beispielsweise scharfen Fal- 
tungskern und einem beispielsweise schmalen Schichtemp- 

^) f indlichkeitsprof il, 

c) Bereitstellen einer Transf erfunktion als funktionaler 
Zusammenhang zwischen Fensterweite und Bildscharfe, 

d) Automatische Berechnung der Bildscharfe des CT-Bildes 
einer ausgewahlten, in dem Primardatensat z liegenden 

30 Schicht in Abhangigkeit von einer ausgewahlten Fensterwei- 

te fur die ausgewahlte Schicht durch einen Bildverarbei- 
tungsprozess auf Basis der Transf erf unktion . 

Dabei wird gemaft einem ersten Ausf uhrungsbeispiel dieses 
35 ersten erf indungsgemafien Verfahrens in dem Bildverarbeitungs- 
prozess die Fensterweite an zumindest einen Parameter des 
dreidimensionalen Faltungskerns gekoppelt, mit dem der Pri- 
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mardatensatz erneut gefaltet und die zumindest eine Schicht 
hinsichtlich der Bildscharfe modifiziert wird. 

Gemafi einen zweiten Ausf uhrungsbeispiel des ersten erfin- 
dungsgemaJien Verfahrens wird in dem Bildverarbeitungsprozess 
die Fensterweite an zumindest einen Parameter eines zweidi- 
mensionalen Faltungskerns gekoppelt, mit dem die zumindest 
eine Schicht gefaltet und hinsichtlich der Bildscharfe modi- 
fiziert wird. 

Bei beiden Ausf iihrungsbeispielen wird eine erneute Rekon- 
struktion mit geanderten Kernparametern und die erneute 
Berechnung von Sekundarschnitten vermieden, was als solches 
ein wesentlich schnelleres Verfahren darstellt. 

Vorteilhaf terweise wird anfangs der Primardatensatz mit einem 
maximal scharfen Faltungskern und maximal schmalen Schicht- 
empf indlichkeitsprof il rekonstruiert . 

20 Erf indungsgemafi kann die Schicht eine axiale Schicht oder 
aber auch eine Sekundarschicht sein. 

Die Auswahl der Fensterweite erfolgt dabei durch den Anwender 
mittels Maus oder Tastatur. 

Weiterhin kann es von Vorteil sein, wenn die Trans f erf unktion 
zusatzlich einen funktionalen Zusammenhang zwischen Fenster- 
weite und Schichtempf indlichkeitsprof il darstellt. 

30 Er f indungsgemaft wird ferner ein zweites erf indungsgemafies 

Verfahren zur Filterung von CT-Bildern vorgeschlagen aufwei- 
send die folgenden Schritte: 



10 



(1) Aufnehmen eines CT-Rohdatensatzes mit einem CT-Gerat oder 
35 mit einem C-Bogen-Gerat , 



200201333 



(2) Rekonstruktion eines Primardatensatzes mit beispielsweise 
scharfem Faltungskern und beispielsweise schmalem 
Schichtempf indlichkeitsprof il, 

(3) Rekonstruktion eines Bildstapels mit einer entsprechenden 
5 Bildcharakteristik auf der Basis des Primardatensatzes, 

(4) Berechnung einer geanderten Bildcharakteristik des Bild- 
stapels durch einen im Hintergrund auf dem Bildrechner 
ablauf enden Bildverarbeitungsprozess, 

(5) Visualisieren des Bildstapels in Form von CT-Bildern mit 
10 der geanderten Bildcharakteristik. 

Dabei kann der Anwender in einem wiederum ersten Ausfuhrungs- 
beispiel dieses zweiten erf indungsgemaften Verfahrens iiber 
eine Eingabeschnittstelle zumindest einen Parameter des 
15 dreidimensionalen Primardaten-Faltungskerns variieren, mit 

dem dann in dem Bildverarbeitungsprozess der Primardatensatz 
erneut gefaltet und daraus ein neuer Bildstapel ermittelt 
wird, der hinsichtlich der Bildscharfe modifiziert wird. 

20 In einem zweiten Ausf uhrungsbeispiel des zweiten erfindungs- 
gemaiien Verfahrens kann der Anwender liber eine Eingabe- 
schnittstelle zumindest einen Parameter eines zweidimensiona- 
len Faltungskerns variieren, mit dem dann in dem Bildverar- 
beitungsprozess die einzelnen Schichten des Bildstapels ein- 

^5 zeln gefaltet und hinsichtlich ihrer Bildscharfe modifiziert 
werden. 

Die Schritte (4) und (5) konnen vom Anwender so oft wieder- 
holt werden, bis eine bef riedigende Bildcharakteristik er- 
30 reicht ist. 

Der Bildstapel kann er f indungsgemaft einen Stapel aus axialen 
Bildern oder aber einen Stapel aus beliebigen Sekundarschnit- 
ten darstellen. 



35 



Ein weiterer Aspekt der vorliegenden Erfindung ist es, dass 
bei ungeniigender Bildcharakteristik auch die Schichtdicke 
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einer in dem Bildstapel enthaltenen Schicht variiert werden 
kann . 



Erf indungsgemafl wird ferner ein Computer-Tomographiegerat zur 
5 Durchfuhrung des Verfahrens nach einem der vorhergehenden 

Anspriiche vorgeschlagen, mit einem Computer, der die einzel- 
nen Schritte der CT-Bildverarbeitung durchfiihrt. 

Ebenso wird er f indungsgemafi ein Computersof tware-Produkt 
10 vorgeschlagen, welches ein Verfahren nach einem der vorherge- 
henden Anspriiche implementiert , wenn es auf einer mit einem 
Computer-Tomographiegerat verbundenen Recheneinrichtung 
lauf t • 



15 Weitere Merkmale, Eigenschaf ten und Vorteile der vorliegenden 
Erfindung werden nunmehr anhand von Ausf uhrungsbeispielen und 
unter Bezugnahme auf die begleitenden Figuren der Zeichnungen 
erlautert . 



20 Figur 1 zeigt schematisch eine CT-Apparatur fur ein Facher- 
strahlverf ahren gemaB der vorliegenden Erfindung, 

Figur 2 zeigt eine Hounsf ield-Skala, in der die Hounsfield- 
Einheiten (HE) unterschiedlicher Organe des mensch- 
^5 lichen Korpers angegeben sind, 

Figur 3 zeigt eine Fensterung (Windowing) bei der Darstel- 
lung von CT-Bildern, 

30 Figur 4 zeigt eine erf indungsgemafte Transf erf unktion, 

Figur 5 zeigt die Aquipotentialf lache eines 3D-Faltungskerns 
in Form eines Ellipsoides, das durch eine charakte- 
ristische Ebene geschnitten wird, 



35 



Figur 6 zeigt die durch die charakteristische Ebene gegebene 
Aquipotentiallinie in Form einer Ellipse. 
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In Figur 1 ist schematisch ein Computer-Tomographiegerat fur 
ein Facherstrahlverf ahren dargestellt, das gemafi der vorlie- 
genden Erfindung arbeitet. Bei diesem Gerat rotieren Rontgen- 
5 Rohre 1 und Strahlenempf anger 2 (Detektoren) gemeinsam urn 

eine Drehmitte, die auch Mitte des kreisf ormigen Messfeldes 5 
ist, und in der sich der zu untersuchende Patient 3 auf einer 
Patientenliege 4 bef indet . Ans telle des Computertomographie- 
Gerats konnte auch ein C-Bogen-Gerat zum Einsatz kommen. Um 

10 verschiedene parallele Ebenen des Patienten 3 untersuchen zu 
konnen, kann die Patientenliege entlang der Korperlangsachse 
verschoben werden. Wie man aus der Zeichnung erkennen kann, 
ergeben sich bei CT-Auf nahmen Transversalschnittbilder , also 
Abbildungen von Korperschichten, die im Wesentlichen senk- 

15 recht zur Korperachse orientiert sind. Diese Schichtdars tel- 
lungsmethode stellt die Verteilung des Schwachungswertes 
M-z(x,y) selbst dar (z ist die Position auf der Korperlangs- 
achse) . Die Computer-Tomographie (im folgenden CT genannt) 
benotigt Projektionen unter sehr vielen Winkeln a. Zur Erzeu- 

20 gung einer Schichtauf nahme wird der von der Rontgenrohre 1 
emittierte Strahlenkegel so ausgeblendet , dass ein ebener 
Strahlenf acher entsteht, der eindimensionale Zentralpro j ekti- 
onen der durchstrahlten Schicht entwirft. Zur exakten Rekon- 
struktion der Verteilung der Schwachungswerte |4. z (x,y) muss 
dieser Strahlenf acher senkrecht auf der Drehachse stehen und 
auiierdem so weit gespreizt sein, dass er aus jeder Projekti- 
onsrichtung a die anvisierte Schicht des Messobjektes voll- 
standig uberdeckt. Dieser das Objekt durchdringende Strahlen- 
facher wird von Detektoren, die auf einem Kreissegment linear 

30 angeordnet sind, aufgefangen. Bei handelsiiblichen Geraten 
sind dies bis zu 1000 Detektoren. Der einzelne Detektor 
reagiert auf die eintref f enden Strahlen mit elektrischen 
Signalen, deren Amplitude proportional zur Intensitat dieser 
Strahlen ist. 

35 

Jedes einzelne zu einer Projektion a gehorige Detektorsignal 
wird jeweils von einer Messelektronik 7 aufgenommen und an 
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einen Computer 8 weitergeleitet . Mit dem Computer 8 lassen 
sich die gemessenen Daten nun in geeigneter Weise verarbeiten 
und zunachst in Form eines Sinugramms (in dem die Projektion 
a als Funktion der Messwerte des entsprechenden Kanals (3 
5 aufgetragen wird) in sogenannten Gordon-Einheiten, schlieft- 
lich aber in Form eines natiirlichen Rontgenbildes in Houns- 
f ield-Einheiten an einem Monitor 6 visualisieren . 

Ziel der vorliegenden Erfindung ist es, einfache und schnelle 
10 Verfahren bereit zustellen, urn die Bildcharakteristik von CT- 
Bildern - insbesondere von Sekundarschnitten - nach diagnos- 
^ tischen Gesichtspunkten anpassen zu konnen . 

Ein erstes Ausf llhrungsbeispiel eines ersten erf indungsgemafien 
15 Verfahrens besteht darin, - im Gegensatz zum oben beschriebe- 
nen Stand der Technik - eine direkte Filterung des aus den 
Rohdaten gewonnenen Volumendatensatzes vorzunehmen. Dabei 
werden die Rohdaten in einem ersten Schritt mittels Detektor- 
elementen 2 in Form von Schwachungsprof ilen unterschiedlicher 
20 Projektionen und diese wiederum in unterschiedlichen benach- 
barten Schichten entlang der Korperachse gemessen. Die Rohda- 
ten werden von der Messelektronik 7 aufgenommen und an den 
Computer 8 weitergeleitet, der in einem zweiten Schritt aus 
diesem CT-Rohdatensat z einen Primardatensatz (Volumendaten- 
satz) errechnet, welcher als solcher in Form von HU-Einheiten 
die Grundlage fiir die Berechnung axialer bzw. sekundarer 
Schnittbilder darstellt. Alle weiteren Schritte des erfin- 
dungsgemaften Verfahrens werden in dem Computer 8 des CT- 
Gerates rein rechnerisch durchgefuhrt bzw. auf dem Monitor 6 
30 visualisiert . 

So erfolgt die erf indungsgemaBe Filterung des Primardatensat- 
zes retrospektiv in dem ersten Ausf uhrungsbeispiel des ersten 
erf indungsgemaiien Verfahrens durch eine dreidimensionale 
35 Faltung, indem ein dreidimensionaler Faltungskern in Form 

einer dreidimensionalen Matrix direkt auf den Volumendaten- 
satz angewendet wird. Der Volumendatensatz wird so in der 
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Weise verandert, dass die daraus erzeugten Sekundarschnitte 
die gewunschte Bildcharakteristik aufweisen. 

Das Frequenzverhalten der axialen Primarbilder in transversa- 
5 ler Richtung (also in x-y-Richtung) ist durch die Modulati- 
ons-Ubertragungs-Funktion des verwendeten Rekonstruktionsal- 
gorithmus bestimmt und im allgemeinen isotrop. Die Modula- 
tions-Ubertragungs-Funktion gibt an, welche (Orts-) Frequenz 
mit welcher Amplitude im Bild sichtbar ist. Dabei ist die 

10 Modulations-Ubertragungs-Funktion im Wesentlichen von der 

Geometrie des Rontgensystems ( Fokusbreite, Anzahl der Detek- 
torkanale in transversaler Richtung, Anzahl der Projektionen 
pro Umlauf usw. ) und von dem Faltungskern, der bei der Rekon- 
struktion verwendet wird, abhangig. In z-Richtung definiert 

15 das SSP bestehend aus Messsystemsanteilen (Fokuslange, Breite 
der Detektorelemente usw.) und Algorithmusanteilen (axiale 
Gewichtungsf unktion im Spiralalgorithmus ) das Verhalten im 
Frequenzraum. In der Regel ist diese Richtung - auch bei 3D- 
isotroper Auflosung - gegenuber den axialen Schichtebenen 

20 ausgezeichnet . 

Ein dreidimensionaler Faltungskern mit der gewiinschten Wir- 
kung auf die Bildcharakteristik kann radialsymmetrisch in 
Bezug auf die z-Achse sein, wahrend die z-Achse eine Sonder- 
^ rolle einnimmt. Im Allgemeinen ist die Filtercharakteristik 
jedoch entlang aller drei Raumachsen unterschiedlich . 

Zur Veranschaulichung sei angenommen, dass der zur Filterung 
korrespondierende 3D-Fal tungskern eine Gauiif unktion sei, wo- 

30 bei im allgemeinen Fall deren Breite entlang aller Achsen 

verschieden ist. Die Aquipotential f lache eines solchen Fal- 
tungskerns ist dann die Oberflache eines Ellipsoides mit 
unterschiedlichen Halbachsen, wie es beispielsweise in Figur 
5 dargestellt ist. Auf Basis eines solchen Ellipsoides wird 

35 der Volumendatensat z ermittelt auf dessen Grundlage dann die 
Sekundarschnitte (MPRs) interpoliert werden. 
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Ein zweites Ausf uhrungsbeispiel des ersten erf indungsgemaiien 
Verfahrens besteht darin, die Sekundarschnitte (MPRs) direkt 
zu filtern. Das ubliche Vorgehen dabei ist, aufgrund der zu- 
vor eingestellten Schichtdicke sowie der raumlichen Lage 
5 eines Sekundarschnittes eine zweidimensionale Faltung vorzu- 
nehmen. Der dazu benotigte zweidimensionale Faltungskern ist 
eine zweidimensionale Matrix, die durch den Schnitt des 3D- 
Faltungskerns mit der charakteristischen Ebene des entspre- 
chenden Sekundarschnittes gegeben ist. Die charakteristische 

10 Ebene ist durch eine zu den Sekundarschnitten parallele Ebene 
durch den Ursprung eindeutig festgelegt und definiert die 

- Lage der Sekundarschnitte. Im allgemeinen Fall ergeben sich 
fur die Aquipotentiallinien des zweidimensionalen Faltungs- 
kernes Ellipsen mit im Allgemeinen gegenuber der kartesischen 

15 Geometrie gedrehten Achsen. In Figur 6 ist eine solche Ellip- 
se dargestellt. 

Im Grunde genommen sind beide Ausf iihrungsbeispiele - 3D- 
Faltung des Volumendatensat zes und anschliefiende Berechnung 

20 der MPRs (3D-Filter + MPR) bzw. Berechnung der MPRs und deren 
anschlieiiende Filterung mit einem 2D-Faltungskern (MPR + 2D- 
Filter) - Equivalent, sofern die Eigenschaf ten senkrecht zur 
MPR-Ebene durch den Algorithmus zur Berechnung der Sekundar- 
schnitte bestimmt werden. So wird beispielsweise bei der Be- 

^ rechnung einer MPR mit Hilfe einer Interpolationsfunktion vom 
Benutzer die Schichtdicke eingestellt. Insbesondere bei einer 
zweidimensionalen Filterung sollen auch nur die transversalen 
Charakteristika der MPRs verandert werden. 

30 Speziell das zweite Ausf uhrungsbeispiel des ersten erfin- 

dungsgemaften Verfahrens ( zweidimensionaler Faltungskern ange- 
wendet auf bereits vorhandene MPRs) ermoglicht die Verande- 
rung der Bildcharakteristik (Scharfe, Rauschen) bereits vor- 
handener MPRs in einem einzigen Schritt. Bei typischen Re- 

35 chenzeiten von ca. lOOmsek pro MPR ist dieses Verfahren we- 
sentlich schneller als das herkommliche Verfahren, bei dem 
aus den Rohdaten eine erneute Rekonstruktion des Primardaten- 
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satzes und eine anschliefiende erneute Berechnung der Sekun- 
darschnitte durchgefiihrt wird. Legt man einen Stapel von 100 
Primarbildern zugrunde, aus dem 50 Sekundarschnitte berechnet 
wurden, so erfordert das konventionelle Verfahren allein fur 
5 die Rekonstruktion 50 Sekunden wenn man zwei rekonstruierte 
Primarbilder pro Sekunde annimmt . Hinzu kommt aufierdem noch 
die Erzeugung der Sekundarschnitte. Eine direkte Filterung 
gemafi des zweiten Ausf iihrungsbeispieles des ersten Verfahrens 
der vorliegenden Erfindung hingegen benotigt nur etwa 5 
10 Sekunden. 

Das Ausfuhrungsbeispiel des ersten erf indungsgemaften Verfah- 
rens in Form der 3D-Faltung ist nicht ganz so schnell wie die 
direkte Faltung der bereits berechneten Sekundarschnitte, da 

15 vor einer erneuten Berechnung von Sekundarschnitten der Volu- 
mendatensatz mit einem 3D-Faltungskern neu bestimmt wird, was 
zeitaufwendiger ist als eine reine Faltung bestehender MPRs . 
Der Vorteil einer solchen Vorgehensweise besteht jedoch da- 
rin, dass im Nachhinein beliebige Sekundarschnitte durch das 

20 neu berechnete Volumen gelegt werden konnen, die letztendlich 
die gewiinschte Bildcharakteristik aufweisen. 

Ein weiterer Aspekt der vorliegenden Erfindung ist es, in 
einem zweiten er f indungsgemaften Verfahren dem Anwender eine 
/jj^5 Moglichkeit zu geben, mittels einer geeigneten visuellen 

Schnittstelle (engl.: Interface), die letztlich im Computer 8 
implementiert werden soil, auf einfache und schnelle Weise 
die Bildcharakteristik eines CT-Bildes ( insbesondere einer 
MPR) der diagnostischen Anforderung anpassen zu konnen. Spe- 
30 ziell konnen so unterschiedliche anatomische Gewebstypen mit 
unterschiedlicher Bildcharakteristik (Bildscharfe und Bild- 
rauschen) dargestellt werden. 

Das zweite erf indungsgemafie Verfahren der vorliegenden Erfin- 
35 dung macht sich die Tatsache zunutze, dass die CT-Werte (HU- 
Werte) je nach der organischen Struktur in der Houndsfield- 
Skala unterschiedliche Bereiche belegen. 
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losen, wahlt man geeigneter Weise ein kleines Fenster, in dem 
die Bildscharfe zwar hoch ist, gleichzeitig aber eine damit 
verbundene hohe Rauschamplitude in Kauf genommen wird. 

5 Aus diagnostischer Sicht ist es sinnvoll, wenn die zur Bild- 
darstellung gewahlte Fensterweite die Schar f eparameter und 
damit die Bildcharakteristik festlegt. Der erf inderische 
Gedanke besteht also darin, die zur Bilddarstellung gewahlte 
Fensterweite an Bildscharfe und Rauschen zu koppeln. Dazu 

10 wird erf indungsgemaft eine Transfer funktion definiert, die der 
gewahlten Fensterweite eindeutig einen Faltungskern zuordnet. 
Eine solche Transf erf unktion ist in Figur 4 dargestellt. Die 

^ Abszisse stellt die Fensterweite, die Ordinate die Bildschar- 
fe dar. Sinnvollerweise ist die Transf erf unktion eine monoton 

15 ansteigende Kurve, da mit zunehmender Fensterweite bei gege- 
bener maximaler Rausch-Amplitude das Rauschen immer weiter 
unterdruckt und der Bildkontrast erhoht wird. Generell jedoch 
kann ein beliebiger, die Diagnose sinnvoll unterstut zender 
Zusammenhang zwischen Fensterweite und Bildscharfe gewahlt 

2 0 werden. 



Durch diesen funktionalen Zusammenhang zwischen Fensterweite 
und Bildscharfe, der im Computer 8 des CT-Gerates erfindungs- 
gemaiJ sof twaremaiiig implementiert ist, ist es dem Anwender 

^5 moglich, gemaft einem ersten Ausf uhrungsbeispiel des zweiten 
erf indungsgemaften Verfahrens die Bildscharfe des CT-Bildes 
einer ausgewahlten, in dem Volumendatensatz liegenden Schicht 
durch Filtern dieses zugrundeliegenden Volumendatensat zes in 
einem Bildverarbeitungsprozess auf Basis der Transf erf unktion 

30 zu andern. Im allgemeinen Fall kann dabei - wie bereits oben 
ausftihrlich beschrieben - die Filtercharakteristik (die Be- 
schaffenheit des der Filterung zugrundeliegenden Faltungs- 
kerns) entlang aller drei Raumachsen unterschiedlich sein. 
Anschaulich sind dann die Aquipotentialf lachen des Faltungs- 

35 kerns Oberflachen von Ellipsoiden mit unterschiedlichen Halb- 
achsen. Ein Filtern des Volumendatensat zes hat den Vorteil, 
dass im Nachhinein beliebige Schichten in diesem Volumenda- 
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tensatz mit der entsprechend modif izierten Bildcharakteristik 
erhalten werden konnen. 

Aus Grlinden der (Rechen-) Geschwindigkeit wird man jedoch in 
5 einem zweiten Ausfuhrungsbeispiel des zweiten er f indungsgema- 
lien Verfahrens eine direkte Filterung der Sekundarschnitte 
durch Verwendung eines zweidimensionalen Faltungskerns vor- 
nehmen. Die raumliche Lage der Sekundarschnitte (MPRs) ist 
durch eine zu den Sekundarschnitten parallelen Ebene durch 
10 den Koordinatenursprung (im Folgenden als "charakteris tische 
Ebene" bezeichnet) eindeutig festgelegt. Der verwendete 2D- 
Faltungskern ist - wie bereits erlautert und in Figur 6 dar- 
gestellt - als Schnitt des 3D- Faltungskerns (Ellipsoid) mit 
der charakteristischen Ebene in Form einer Ellipse gegeben. 

15 

Die Filterung nach dem zweiten Ausfuhrungsbeispiel des zwei- 
ten erf indungsgemaBen Verfahrens wird dadurch realisiert, 
dass eine zweidimensionale Faltung auf der entsprechenden 
Sekundarebene (MPR) eines vorgegebenen Ref erenzvolumens 
20 durchgefuhrt wird. Das Ref erenzvolumen ist der Primardaten- 
satz (ursprunglicher Volumendatensatz ) , der beispielsweise 
mit maximal scharfem Faltungskern sowie maximal schmalem 
Schichtempf indlichkeitsprof il aus dem Rohdatensatz errechnet 
wurde . 

Allgemein wird eine obige Filterung nach dem ersten Ausfuh- 
rungsbeispiel des zweiten erf indungsgemafien Verfahrens so 
durchgefuhrt, dass nach Vorgabe der Transf erf unktion die 
dreidimensionale Filtercharakteristik in Abhangigkeit von der 

30 Fensterweite durch einen Bildverarbeitungsprozess , der auf 
dem Computer 8 durchgefuhrt wird, adaptiert wird. Das Refe- 
renzvolumen wird mit dem entsprechenden dreidimensionalen 
Faltungskern der gewunschten Bildcharakteristik angepasst und 
daraus anschlieftend die entsprechenden Sekundarebenen (MPRs) 

35 berechnet. Schlieftlich wird die jeweilige Sekundarebene mit 
der gewahlten Fensterweite zur Darstellung gebracht . 
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Die praktische Realisierung einer derartigen fenstergesteuer- 
ten Filterung von CT-Bildern erfolgt beispielsweise dadurch, 
dass der Anwender auf dem Monitor 6 des Computers 8 auf einem 
virtuellen Interface mit der Maus die Fensterweite (Fenster- 
5 einstellung) entsprechend den diagnostischen Zwecken andert 
und gemafi der Transf erf unktion die Filtercharakteristik des 
zugrundeliegenden Faltungskerns entsprechend zeitgleich ver- 
andert. So lasst sich allein durch Variation der Fensterweite 
die Bildcharakteristik (speziell Scharfe und Rauschen) den 
10 diagnostischen Bediirfnissen des Anwenders anpassen. 

Beide beschriebenen erf indungsgemafien Verf ahren sind nicht 
nur fur die adaptierte Darstellung von Sekundarschnittbildern 
sondern auch in der klinischen Routine im Besonderen fur die 
15 diagnostische Auswertung von axialen Schichtbildern bedeut- 
sam. 
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Patentanspruche 

1. Verfahren zur f enstergesteuerten Filterung von CT-Bildern 
aufweisend die folgenden Schritte: 

5 

a) Aufnehmen eines CT-Rohdatensatzes mit einem CT-Gerat oder 
mit einem C-Bogen-Gerat , 

b) Rekonstruktion eines Primardatensat zes aus dem CT- 
Rohdatensatz mittels einem beispielsweise scharfen Fal- 

10 tungskern und einem beispielsweise schmalen Schichtemp- 

f indlichkeitsprof il, 
^ c) Bereitstellen einer Transf erf unktion als funktionaler 
Zusammenhang zwischen Fensterweite und Bildscharfe, 
d) Automatische Berechnung der Bildscharfe des CT-Bildes 
15 einer ausgewahlten, in dem Primardatensatz liegenden 

Schicht in Abhangigkeit von einer ausgewahlten Fensterwei- 
te fur die ausgewahlte Schicht durch einen Bildverarbei- 
tungsprozess auf Basis der Transf erf unktion . 

20 2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, 
dass in dem Bildverarbeitungsprozess die Fensterweite an 
zumindest einen Parameter des dreidimensionalen Fal tungskerns 
gekoppelt ist, mit dem der Primardatensatz erneut gefaltet 

£5 wird und die zumindest eine Schicht hinsichtlich der Bild- 
scharfe modifiziert wird. 

3. Verfahren nach Anspruch 1, 
da durch gekennzeichnet, 
30 dass in dem Bildverarbeitungsprozess die Fensterweite an 

zumindest einen Parameter eines zweidimensionalen Faltungs- 
kerns gekoppelt ist mit dem die zumindest eine Schicht gefal- 
tet und hinsichtlich der Bildscharfe modifiziert wird. 



35 



4. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, 
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10. Verfahren nach Anspruch 9, 
dadurch gekennzeichnet, 

5 dass der Anwender uber eine Eingabeschnittstelle zumindest 
einen Parameter des dreidimensionalen Primardaten- 
Faltungskerns variieren kann, mit dem dann in dem Bildverar- 
beitungsprozess der Primardatensat z erneut gefaltet und 
daraus ein neuer Bildstapel ermittelt wird der hinsichtlich 
10 der Bildscharfe modifiziert ist. 

11. Verfahren nach Anspruch 9 bis 10, 
dadurch gekennzeichnet, 

dass der Anwender liber eine Eingabeschnittstelle zumindest 
15 einen Parameter eines zweidimensionalen Faltungskerns variie- 
ren kann, mit dem dann in dem Bildverarbeitungsprozess die 
einzelnen Schichten des Bildstapels einzeln gefaltet und 
hinsichtlich ihrer Bildscharfe modifiziert werden. 

20 12. Verfahren nach Anspruch 9 bis 11, 

dadurch gekennzeichnet, 

dass die Schritte (4) und (5) vom Anwender so oft wiederholt 
werden konnen, bis eine bef riedigende Bildcharakteristik 
erreicht ist. 

13. Verfahren nach Anspruch 9 bis 12, 
dadurch gekennzeichnet, 

dass der Bildstapel einen Stapel aus axialen Bildern dar- 
stellt . 

30 

14. Verfahren nach Anspruch 9 bis 12, 
dadurch gekennzeichnet, 

dass der Bildstapel einen Stapel aus beliebigen Sekundar- 
schnitten darstellt . 

35 

15. Verfahren nach Anspruch 9 bis 14, 
dadurch gekennzeichnet, 
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dass bei ungenugender Bildcharakteristik auch die Schichtdi- 
cke einer in dem Bildstapel enthaltenen Schicht variiert 
wird. 

16. Computer-Tomographiegerat zur Durchfiihrung des Verfahrens 
nach einem der Anspruche 1 bis 15, mit einem Computer, der 
die einzelnen Schritte der CT-Bildverarbeitung durchfuhrt. 

17 . Computersof tware-Produkt , 

dadurch gekennzeichnet, 
dass es ein Verfahren nach einem der Anspruche 1 bis 15 
implementiert , wenn es auf einer mit einem Computer- 
Tomographiegerat verbundenen Recheneinrichtung lauf t . 
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Zusainmen fas sung 

Retrospektive bzw. f enstergesteuerte Filterung von Bildern 
zur Adaption von Scharfe und Rauschen in der Computer- 
Tomographie 

Die vorliegende Erfindung bezieht sich auf ein Verfahren zur 
beispielsweise f enstergesteuerten Filterung von CT-Bildern 
aufweisend die folgenden Schritte: 

a) Aufnehmen eines CT-Rohdatensatzes mit einem CT-Gerat oder 
mit einem C-Bogen-Gerat , 

b) Rekonstruktion eines Primardatensatzes aus dem CT- 
Rohdatensatz mittels einem beispielsweise scharfen Fal- 
tungskern und schmalen Schichtempf indlichkeitsprof il, 

c) Bereitstellen einer Transf erf unktion als funktionaler 
Zusammenhang zwischen Fensterweite und Bildscharfe, 

d) Automatische Berechnung der Bildscharfe des CT-Bildes 
einer ausgewahlten, in dem Primardatensatz liegenden 
Schicht in Abhangigkeit von einer ausgewahlten Fensterwei- 
te fur die ausgewahlte Schicht durch einen Bildverarbei- 
tungsprozess auf Basis der Transf erf unktion . 
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